Journal of Physics A. Vol.38. No.10/11. (2005) pp.2145-2155. 



Phase-Space Metric for Non-Hamiltonian Systems 

Vasily E. Tarasov 

Skobeltsyn Institute of Nuclear Physics, Moscow State University, Moscow 119992, Russia 

E-mail: tarasov ©theory, sinp.msu. ru 



We consider an invariant skew-symmetric phase-space metric for non-Hamiltonian sys- 
tems. We say that the metric is an invariant if the metric tensor field is an integral of 
motion. We derive the time-dependent skew-symmetric phase-space metric that satisfies 
the Jacobi identity. The example of non-Hamiltonian systems with linear friction term is 
considered. 
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1 Introduction 

The dynamics of Hamiltonian systems is characterized by conservation of phase-space volume 
under time evolution. This conservation of the phase volume is a cornerstone of conventional 
statistical mechanics of Hamiltonian systems. At a mathematical level, conservation of phase- 
space volume is considered as a consequence of the existence of an invariant symplectic form 
(skew-symmetric phase-space metric) in the phase-space of Hamiltonian systems (HEUHj. 

The classical statistical mechanics of non-Hamiltonian systems is of strong theoretical in- 
terest |llElinilZllSliniiriIlinilI2llinilIlliniliniin|. Non-Hamiltonian systems have been used 
in molecular dynamics simulation to achieve the calculation of statistical averages in various 
ensemble jSJ EJ 112 UH] , and in the treatment of nonequilibrium steady states [13 H3 EH 120] • 
Non-Hamiltonian systems are characterized by nonzero phase space compressibility, and the 
usual phase-space volume is no longer necessarily conserved. 

Tuckerman et al. have argued E] that there is a measure conservation law that involves a 
nontrivial phase-space metric. This suggests that phase-space should be carefully treated using 
the general rules of the geometry of manifolds PUI2|- Tuckerman et al. have applied the concepts 
of Riemannian geometry to the classical statistical mechanics of non-Hamiltonian systems |3J 
E1E]- Tuckerman et al. have argued that, through introduction of metric determinant factors 



g(x, t), it is possible to define an invariant phase-space measure for non-Hamiltonian systems. 



In their approach the metric determinant factor Jg{pc,t), where g(x.,t) is the determinant of 
the metric tensor, is defined by the compressibility of non-Hamiltonian systems. However 
Tuckerman et al. consider only the determinant g(x, t) of the metric. The phase-space metric 
is not considered in [U El E] • Note that Tuckerman et al. suppose that the metric determinant 
factor is connected with symmetric phase-space metric. It can be proved that the proposal to 
use an invariant time-dependent metric determinant factor in the volume element corresponds 
precisely to finding a skew-symmetric phase-space metric (symplectic form) that is an integral 
of motion. Therefore we must consider the skew-symmetric phase-space metric. 
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Sergi fUJ has considered an antisymmetric phase-space tensor field, whose elements are 
general function of phase-space coordinates. In [HH ITT] , the generalization of Poisson brackets 
for the non-Hamiltonian systems was suggested. However the Jacobi identity is not satisfied 
by the generalized brackets and skew-symmetric phase-space metric. As a result the algebra of 
phase-space functions is not time translation invariant. The generalized brackets do not define 
a Lie algebra in phase-space. Note that the generalized brackets of two constant of motion is 
no longer a constant of motion. 

In the present paper we consider an invariant skew-symmetric (antisymmetric) phase-space 
metric for non-Hamiltonian systems. We say that the metric is an invariant, if the metric 
tensor field is an integral of motion. We define the phase-space metric such that the Jacobi 
identity is satisfied. The suggested skew-symmetric phase-space metric allows us to introduce 
the generalization of the Poisson brackets for non-Hamiltonian systems such that the Jacobi 
identity is satisfied by the generalized Poisson brackets. As a result the algebra of phase-space 
functions is time translation invariant. The generalized Poisson brackets define a Lie algebra in 
phase-space. The suggested Poisson brackets of two constant of motion is a constant of motion. 

In section 2, the definitions of the antisymmetric phase-space metric, mathematical back- 
ground and notations are considered. In section 3, we define the non-Hamiltonian systems, and 
consider the Helmholtz conditions. In section 4, we consider the time evolution of phase-space 
metric. We derive the phase-space metric that is an integral of motion. In section 5, the gen- 
eralized Poisson brackets for non-Hamiltonian systems are defined. In section 6, the example 
of phase-space metric for non-Hamiltonian system with the linear friction term is considered. 
Finally, a short conclusion is given in section 7. 

2 Phase-space metric 

The 2n-dimensional differentiable manifold is denoted M. Coordinates are x = (x 1 , x 2n ). 
We assume the existence of a time-dependent metric tensor field Uki(x.,t) on the manifold M. 
We can define a differential 2-form 

u = u; fc /(x, t)dx h A dx\ (1) 

where oj k \ = u k i(x.,t) is a skew-symmetric tensor u kl = —ui k , and the tensor elements a>jy(x, t) 
are explicit functions of time. Here and later we mean the sum on the repeated index k and I 
from 1 to 2n. 

We suppose that the differential 2-form u is a closed nondegenerated form: 

1) If the metric determinant is not equal to zero 

ff (x,t) = det(a; j u(x,t))^0 (2) 

for all points x G M, then the form u is nondegenerated. 

2) If the Jacobi identity 

d k Lj lm + diLj mk + d m uj kl = 0, d k = d/dx k (3) 
for the metric u k i = Wjw(x, t) is satisfied, then the differential 2-form to is closed (du = 0). 
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Phase-space is therefore assumed to be a symplectic manifold. 



Definition 1. A symplectic manifold is a differentiable manifold M with a closed non- 
degenerated differential 2- form to. 



For symplectic manifold, we have the phase-space volume element 



v = —uj n = —uj A ... Aw, 



This differentiable 2n-form can be represented by 

v = ^Jg(x,t)dx 1 A ... A dx 2n , 

where gr(x, t) is defined by equation (J2J). The nondegenerated condition (<?(x, t) 7^ 0) for the 
metric 0Jki{~ x -, t) is equivalent to the condition uj n ^ or v ^ 0. 

It is known [2, that there exists the local coordinates (q, p) such that 

uj = dijdq' A dp 3 . (4) 

Here and later we mean the sum on the repeated index i and j from 1 to n. 



3 Non-Hamiltonian system 

The dynamics is described by a smooth vector field X = X(x), 

with components X k in basis dk = d/dx k . For simplicity, we consider the case where the vector 
field X is time independent. In local coordinates {x k }, equation (J3J) has the form 

^ = X k . (6) 
dt v ; 

Consider now the definition of the Hamiltonian systems, which is used in jj. 

Definition 2. A classical system on the symplectic manifold (M, uj) is called a Hamil- 
tonian system if the differential 1-form uj{X) is a closed form 

du(X) = 0, 

where u(X) = ix^ is the contraction (interior product) of the 2-form uj with vector X, and d 
is the exterior derivative. 

A classical system (0j on the symplectic manifold (M,uj) is called a non-Hamiltonian 
system if the differential 1-form w(X) is nonclosed du(X) ^ 0. 
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Proposition 1. The classical system J3J) is a Hamiltonian system if the conditions 

J kl (u, x, t) = d k (u lm X m ) - d^X™) = (7) 

are satisfied. 

Proof. In the local coordinates {x k }, we have 

lo(X) = X k dx k = u M X l dx k , 
where X k = uJ k \X l . In this case, the exterior derivative of 1-form uj(X) is 

duj(X) = d(X k dx k ) = d t X k dx l A dx k . 
Using da A db = —db A da, we get 

dcu(X) = i(<« - diX^dx' A dx l . 

As the result we have the differential 2-form 

duj(X) = ^(d k (oo lm X m ) - d l {u Jkm X m ))dx k A dx l . (8) 

This differential 2-form is a symplectic form, which can be called " non-Hamiltonian symplectic 
form". If the Helmholtz conditions (J7J) are satisfied, then the differential 1-form uj(X) is closed 
(duj(X) = 0), and the classical system © is a Hamiltonian system. 

Let us consider the canonical coordinates x = (x 1 , x n , x n+1 , x 2n ) = (q 1 , ...,q n ,p , ...,p n ). 
Equation (jOJ) can be written as 

| = G%P) , # = F%P) . (9 ) 

Corollary. If the right-hand sides of equations (GJ) satisfy the Helmholtz conditions \21\ [EE/ 
for the phase-space with ^j, which have the following form 

0, (10) 

0, (11) 
0, (12) 



dG l 




dpi 


dp i 


dGi 


8F i 


dq { 


dpi 




dFi 


dqi 


dq l 



then classical system (GJ) is a Hamiltonian syst 



em. 



Proof. In the canonical coordinates {q,p), the vector field X has the components (G l ,F l 
which are used in equation (jHJ). The 1-form uj(X) is defined by the following equation: 

u{X) = ±{G i dp i -F i dq i ), 



where Gi = 5ijG^ and Fi = 5ijFK The exterior derivative for this form can now be written by 
the relation 

du(X) = -(diGidp 1 ) - d{F i dq i 
It now follows that 

1 / r)G F)G f)F f)F \ 

dco(X) = - l-^dqi A dp' + -^dpi A dp' - -^dq j A dg* - -^-dp> A dA . 

This equation can be rewritten in an equivalent form 
, , v . 1 (dG j dFA , i , , 1 fdGj dGA , i ,4 1 fdFi dFA ■ . 

Here we use the skew-symmetry of dq 1 A dqi and dp 1 A cfp 7 with respect index % and j. It is 
obvious that conditions (fTU|) - (|T2*j) lead to the equation du(X) = 0. 



4 Time evolution of phase-space metric 

Let us find a time- dependent symplectic 2-form u that satisfies the equation du/dt = 0. 
It is known 0] as the following proposition. 

Proposition 2. If the system x = X on the symplectic manifold (M, u) with time- 
independent symplectic form (du k i/dt = 0) is a Hamiltonian system, then differential 2-form 
to is conserved, i.e., du/dt — 0. 

Proof. The proof of this theorem is considered in [210]. 

Let us consider a generalization of this proposition. 

Proposition 3. If the time- dependent metric uj k \ = u k i(^,t) is a skew- symmetric metric 
(uiki = —uJik) that is satisfied by the Jacobi identity J3J), and the system is defined by equation 
(OJ), then the total time derivative of the differential 2-form (QJ) is given by 

du = ^dm _ g k{uJimXm) + dl ( Ukm X m ))dx k A dx 1 . (13) 

Proof. The time- derivative of the time-dependent symplectic form uj is given by 

du d / . . , u , a dujki , u , , , ( dx k \ , , , k , ( dx l \ 

— = —(uj k i(x, t)dx A dx = ——dx A dx + u M d — — A dx + u M dx A d —— . 
dt dt y 'at \ at J \dt J 



Then, using the equation 



duj H _ du k i dujki dx m _ du k i 

~w-^r + dx^^f-^r +x dMh 
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and equation (JBJ), we find that 



<L ' ' (kCU ' X m d m uj kl ) dx k A dx l + oo kl dX k A dx l + uj M dx k A dX l . 



dt \ dt 
Using dX k = d m X k dx m , we have 
duj ( duiM 



, + X m d m u kl dx k A dx l + u kl d m X k dx m A dx l + u kl d m X l dx k A dx m . 
at \ at ) 

This expression can be rewritten in an equivalent form 

Tt = + XmdmUkl + U ™ ldkXm + "krndlX™^ dx k A dxK 

Using the rule of term-by-term differentiation in the form 

^mid k X m = d k (uj m iX m ) — X m d k ui m i; u) km d[X m = di(u km X m ) — X m diu km , 
we get the following equation: 

li = + ~ - + + d '^ X ">) A iX '- 

Using the Jacobi identity (jHJ), and skew symmetry oj rn i = —uJi m , uj km = —u mk , we have 

^ = (^f ~ 9 MmX m ) + d^X™)^ dx k A dx l . 

As the result, we obtain equation for the total time derivative of symplectic form. 

The total time derivative of the symplectic form is defined by equation (|T3j). If the total 
derivative is zero, than we have the integral of motion or invariant. It is easy to see that the 
differentiable 2-form uo is invariant if the phase-space metric wjy(x, t) is satisfied by the equation 

^ = d k (u lm X m ) - dt(u; km X m ). (14) 
at 

This equation can be rewritten in an equivalent form 

dojhi 

Kl jms 

Q t — J k l ^msi 

where the operator J is defined by the equation 

47 = \{{5?dk - 5Tdi)X s - (6?d k - SiddX"). (15) 

Proposition 4. The differentiable 2-form uj is invariant (is an integral of motion for 
non-Hamiltonian system $)) if the phase-space metric u k i(x.,t) is defined by the equation 

w«(x,t) = [exp(t J) J w ms (x,0). (16) 
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Here J is an operator that is defined by equation 

Proof. Let us consider the formal solution of equation (jl4j) in the form 

oo jn 

^fci(x,t) = -r4?( x )> 

where = —oj\ k \ In this case, the time-independent tensor fields uj^\x) are defined by the 
recursion relation 

(n+l) _ rs ( (n) y m\ o / M ym \ 

This equation can be rewritten in an equivalent form 

= {5Td k -5Td l ){X s ^l). 
Using the skew symmetry of the ujj^J, we have 
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(n+l) 
kl 



\{{5Td k - 8 k ^d l )x s - (5?d k - s^x^l 



2 

This relation can be represented in the form 

(n+l) _ jms (n) 
u kl — J kl ^msi 

where the operator J is defined by equation f)15|). Therefore the invariant phase-space metric 
is defined by the following equation: 

oo ±n 

u; k i(x,t) = —\ jn ) kl ^ms(x, 0) = (exp(t J) J W ms (x, 0). 

n=0 ^" 

As the result we have equation (fTBj) . 

Comments. We prove that equation (|14j) for the elements of phase-space metric can be 
expressed in terms of an operator J such that 

duj ki _ f ms n? x 

—gjT — J kl ^ms, {*-<) 

where J k f is defined by equation (fT3j) . In this case, the time evolution of the phase-space metric 
from initial condition u) k i(x.Q, 0) to a value w k i(x, t) at the time t can be written by the equation 



u k i(-x,t) = (exp(t J))™ S fj ms (x,0). (18) 

Here the matrix exponential operator, exp(tJ), can be called the metric propagator. The 
operator J can be considered as a metric analog of the Liouville operator. Note that we can 
use the canonical coordinates (q,p) for t = and the coordinate-independent initial metric: 
uj k i{x, 0) = ujfi = const, d m uj k i('x, 0) = 0. Introducing time step At = t/N, we get 

exp(t J) J u; ms (x, 0) = [\exp(At J)] J w TOS (x, 0). 



It is natural to approximate the short-time propagator exp(tJ) using the Trotter theorem 
|23| 1^1 125] . The vector field X = X(x) can be represented in the form 

X = Xi + X2, 

where Xi is a Hamiltonian term such that 

and X 2 is a friction (non-Hamiltonian) term. As the result we have 

(£s*)« = liWh ~ S£di)Xy ' (S?d k ~ S' k di)X^). (19) 

An important consequence would be the ability to formulate rigorous numerical integration 
algorithms based on Trotter- type splittings [23 12H EH] of the classical propagator exp(AtJ). 
The metric propagator for the small time step is 

exp[At(J x + J 2 )] = exp[^At J 2 ]exp[AtJ 1 }exp[^At J 2 ] + 0(At 3 ). 
Finally, we obtain 

exp[tJ] = exp[t{Ji + J 2 )} = (exp\^-AtJ 2 }exp[AtJ l }exp\^-AtJ 2 1 \\ N + 0(At 3 ). 

5 Poisson brackets for non-Hamiltonian systems 

Let us consider the skew-symmetric tensor field u kl = a; (x, t) that is defined by the equations 

ui kl (x, t)u lm (x, t) = u lk (x,t)uj ml (x,t) = Sf. 
As the result this tensor field satisfies the Jacoby identity 

It follows from the Jacoby identity for u^i- 

In proposition 3, we suggest the time-dependent phase-space metric ^(x, t), which satisfies 
the Jacoby identity. As the result we have Lie algebra that is defined by the following brackets: 

{A,5} = W H (x 1 i)^ l B. (20) 

It is easy to prove that these brackets are Poisson brackets. 

In the general case, the rule of term-by-term differentiation with respect to time that has 
the form 

j t {A,B} = {A,B} + {A,B}, (21) 
where A = dA/dt, is not valid for non-Hamiltonian systems. In the general case, we have 

j t {A,B} = {A,B} + {A,B} + J(A,B) 
8 



where 

J{A,B) =u% ) {x,t)d k Ad l B, u$(x,t) =u km J s {d s X m -d m X s ). 

Note that time evolution of the Poisson brackets (}2T)j) for non-Hamiltonian systems can be 
considered as t-deformation [2H] of the Lie algebra in phase-space. 

If we use the invariant phase-space metric, then the rule (}2~T|) is valid. As the result the 
suggested Poisson brackets (j2*Uj) of two constants of motion is a constant of motion and rule 
(|2*T|) is satisfied. 



6 Example: System with linear friction 

Let us consider the non-Hamiltonian system that is described by the equations: 

dq i dH dp 1 dH 



dt dp 1 ' dt dq l 



K)(t)p>, (22) 



where H = T(p) + U(q) and i; j = 1, ...,n. Here T(p) is a kinetic energy, U(q) is a potential 
energy. The term —Kj(t)pi is a friction term. Usually, this system is described by the phase- 
space metric Uki that has the form 

I \uki\ I = ^ _qt q J j (23) 

where the matrix G = \\g-ij\\ is equal to indentity matrix E = Here G T is transpose 

matrix for the matrix G. The symplectic form is defined by equation (J3J) in the form 



uj = 2gijdq l A dpi = b^dq 1 A dp 
where = (1/2)5^-. The phase-space compressibility 



Qqi Qpi 

of the system (}2*2"|) is defined by the spur of the matrix K — \\K^\\ in the form 

K=-Sp\\K i 1 \\. 



Proposition 5. The invariant phase-space metric for the classical system \2ty) has the form 

G(t 
-G T (t) 



Mt)\\=( X* G{ * ] ), (25) 



where the matrix G is defined by 

rt 
'to 



G(t) = G(t )exp f || K]{t) 1 1 dr. (26) 

J tn 
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Proof. Suppose that phase-space metric depends on time t. Therefore the matrix G and 
elements g^i are the functions of the variable t: G — G(t), gij = Qijit). Let us consider the total 
time derivative of the symplectic form 

lu = 2g ij (t)dq l A dp? . 

As the result, we have 

In order to have the invariant phase-space metric (duj/dt = 0), we use the following equation: 

g im K. 



dgij 



dt 

The matrix G = \ \gij\ \ is satisfied the matrix equation 

dG(t) 



dt 



G(t)K(t), (27) 



where K(t) = \ \Kj(t)\ \ is a matrix of friction coefficients. Suppose that G(t ) = E, where E is 
the identity matrix. The solution of equation (J27j) has the form 



G{t) = G{t )exp [ K{r)dr = exp f K{r)dr. (28) 

J to J to 



If the matrix K is diagonal matrix with elements Kj(t) = Kj(t)8p then we have the matrix 
elements t 

gij(t) = dijexp / Kj(r)dr. 



t 



As the result, the invariant phase-space metric for system (J22j) is defined by equations (J25|) and 
The metric determinant factor [3J E] 

\Jg(x,t) = yjdet\\u kl (x,t)\\ 



for the phase-space metric is defined by the relation ^Jg = \J det{GG T ) . Using det (AB) 
(det A) (det B) and detA T = det A, we have the metric determinant factor in the form ^fg 
\detG\. If the matrix G has the form (|28jl . then 



g[t) = \det exp [ K(r)dr 

Jtn 



Using det exp A = exp Sp A, we get that the invariant phase-space metric has the determinant 
that is connected with the phase-space compressibility jU E] by the equation 

g(i) = exp / SpK(r)dr = exp — / K(r)dr, (29) 

J to Jto 

where k is the phase-space compressibility 



10 



For example, the system 



dqi 
dt 

dq 2 
dt 



Pi 

m ' 

P2 

m ' 



dpi 
dt 

dp 2 
dt 



dU(q) 
dqi 

dU(q) 
dq 2 



- K2P2, 



has the invariant phase-space metric Uki(t) in the form 

\\uki(t)\\ 



1 








e K lt 
















e K 2 t 




_ e K lt 











\ 





_ e K 2 t 









(30) 



In this case, the metric determinant factor is equal to 



yjgify = y/det\\u k i\ 



7 Conclusion 

Tuckerman et al. j3J EJ |B] suggest a formulation of non-Hamiltonian statistical mechanics which 
uses the invariant phase-space measure. The invariant measure is connected with phase-space 
metric. Tuckerman et al. consider the properties of only the determinant g(x, t) of the metric. 
The phase-space metric is not considered in jU E] . We consider the invariant phase-space 
metric for non-Hamiltonian systems. The proposal to use an invariant time-dependent metric 
determinant factor in the volume element corresponds precisely to finding a skew-symmetric 
phase-space metric (symplectic form) that is an integral of motion. Therefore we consider the 
skew-symmetric metric. 

Sergi pm ITT] uses the skew-symmetric phase-space metric that is not satisfied by the Ja- 
coby identity. As the result the generalization of the Poisson brackets for non-Hamiltonian 
systems leads one to non-Lie algebra. Note that non-Lie algebra for non-Hamiltonian sys- 
tems is considered in In this paper we consider an invariant antisymmetric phase-space 
metric that satisfies the Jacoby identity, and defines the Lie algebra in phase-space. We 
call the metric is invariant if the metric tensor field is an integral of motion (du/dt = 0). 
This invariant phase-space metric w&j(x, t) defines the invariant phase-space measure v by the 
equation v = (l/n\)u n = \Jg(x., i)dx x A ... A dx 2n , where g(x,t) is the metric determinant 
g(x, t) = det(u)ki(x.,t)) and dv/dt = 0. The suggested time-dependent skew-symmetric phase- 
space metric leads to a constant value of the entropy density, so that the associated distribution 
function obeys an evolution equation associated with incompressible dynamical flow. 

Note that the invariant phase-space metric of some non-Hamiltonian systems can leads us 
to the lack of smoothness of the metric. In this case, the phase-space probability distribution 
can be collapsed onto a fractal set of dimensionality lower than in the Hamiltonian case [251 l2Zj • 
Unfortunately the description of lack of smoothness in 12Z| is considered without using the 
curved phase-space approach jU • Note that classical systems that are Hamiltonian systems 
in the usual phase-space are non-Hamiltonian systems in the fractional phase-space [T^l 13*2]. 
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In the papers [3U1 EH] , the quantization of the evolution equations for non-Hamiltonian and 
dissipative systems was suggested. Using this quantization it is easy to derive the quantum 
analog of the invariant Poisson brackets, which satisfy the rule of term-by-term differentiation 
with respect to time. 
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